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Abstract 

A nonuniform Neumann boundary-value problem is considered for the Poisson equa¬ 
tion in a thin domain coinciding with two thin rectangles connected through a joint 
of diameter 0(e) . A rigorous procedure is developed to construct the complete asymp¬ 
totic expansion for the solution as the small parameter e —^ 0. Energetic and uniform 
pointwise estimates for the difference between the solution of the starting problem 
(e > 0) and the solution of the corresponding limit problem (e = 0) are proved, from 
which the influence of the geometric irregularity of the joint is observed. 
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Introduction 

Investigations of various physical and biological processes in channels are urgent for numerous 
fields of natural sciences. Special interest of researchers is focused on various effects observed 
in vicinities of local irregularities of the geometry (widening or narrowing) of channels (e.g., 
adhesion to the walls, welds, and stenosis). Results of recent theoretical, experimental and 
numerical studies of flows and wall-pressure fluctuations in channels with different types of 
narrowing are summarized in [H |2l |3] and references therein. Also the study of influence of 
local geometrical irregularities is very important in engineering, since such irregularities often 
directly affect the strength (stability, resistance, power, etc.) of constructions and devices. 
A fairly complete review on this topic has been presented in the article by Gaudiello and 
Kolpakov [1]. 

In addition, the paper [4] is a pioneering paper, where the influence of a local geometrical 
irregularity in a thin domain was studied with the help of multi-scale approach. Videlicet, 
the authors derived the limit problem for a homogeneous Neumann problem for the Poisson 
equation with a right-hand side that depends only of one longitudinal variable in junctions of 
thin domains and showed that the local geometric irregularity in the joint does not affect the 
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view of the corresponding limit problem. However, the convergence theorem and asymptotic 
estimates have not been proved. 

It should be stressed that the error estimates and convergence rate are very important 
both for the justihcation of the adequacy of one- or two-dimensional models aimed at the 
description of actual three-dimensional thin bodies and for the study of boundary effects and 
effects of local (internal) inhomogeneities in applied problems. Particular importance for 
engineering practice is pointwise estimates for approximations, since large values of tearing 
stresses in small region at hrst involve local material damage and then the destruction of 
whole construction. Those estimates can be obtained and substantiated as a result of the 
development of new asymptotic methods. 

In [5j we proved the error estimates and constructed the asymptotic expansion for the 
solution of a boundary-value problem in a thin cascade domain without joints (see Fig. [1]). 
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Figure 1: Thin cascade domain without local joints 


The present paper is devoted to further development of the asymptotic method proposed 
in [5]. Namely, we consider a nonuniform Neumann boundary-value problem for the Poisson 
equation with the right-hand side that depends both on longitudinal and transversal variables 
in a thin cascade domain that consists of two thin rectangles of different thicknesses and 
local geometric irregularity (the joint) between them (this can be either a local widening 
(see. Fig. |2]) or local narrowing (see Fig. |3])). As we will see, there is no essential difference 
between the construction of the asymptotic expansion for the solution to a boundary-value 
problem in 2 - or 3 - dimensional thin cascade domain formed by two thin rods. Therefore, 
to simplify calculations, we consider the two-dimensional case. 



A principal new feature of this paper in comparison with the paper [1] is the construction 
and justihcation of the complete asymptotic expansion for the solution and the proof both 
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Figure 3: Thin cascade domain with a local narrowing 


energetic and pointwise uniform estimates for the difference between the solution of the 
starting problem [e > 0) and the solution of the corresponding limit problem (e = 0). As a 
result, it became possible to identify more precisely the impact of the geometric irregularity 
and material characteristics of the joint on some properties of the whole structure. In 
addition, on the one hand, our results conhrm and complement some of conclusions of the 
article [1], on the other hand show that the main assumptions made in this article are not 
correct. A more informative discussion is given in the last section of the present paper. 

The paper is organized as follows. In Section 2, we construct the formal asymptotic 
expansion for the solution to the problem ([T]). To perform this we generalize the asymptotic 
method for boundary-value problems in thin domains with constant thickness proposed in 
monograph [9]. In particular, we introduced a special inner asymptotic expansion in a 
neighborhood of the joint, determine its coefficients and study some their properties as 
solutions to corresponding boundary-value problems in an unbounded domain. Thus, the 
asymptotics for the solution consists of three parts: the regular part, the boundary parts 
near the extreme vertical sides and the inner part in a neighborhood of the joint. 

In Section 3, we justify the asymptotics fTheorem 13.ip and prove asymptotic estimates 
for the leading terms of the asymptotics fCorollary 13.ip . In Conclusions we analyze results 
obtained in this paper and discuss possible generalizations. 

1 Statement of the problem 

The model thin cascade domain consists of two thin rectangles 

Sjy = ((-1, -ft) X yl) and SlfI = ((at 1) x Tf>) 

that are joined through 12^°^ (referred in the sequel ’’joint”). Here , * = 

1,2; £ is a small parameter; /, hi and h 2 are hxed positive constants. 

The joint are formed by the homothetic transformation with the coefficient £ from 
a domain , i.e., In addition, we assume that 

bl < £ max{hi,h 2 }} C {(a;,|/) : \x\ < s^-)} 

and the interior of the union _ _ _ 

U U 

is a domain with the Lipschitz boundary, which we denote by (see e.g. Fig. 0]). 
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Figure 4: The model thin cascade domain hlg with a local geometric irregularity 


Remark 1.1. As an example of the joint we can consider the following domain: 

^(0) = ^^(0) ^ [{x,y) : -ehg (f) <y < ehf{f), y e (-e^, e^)}, 

where : —h-o (0 <V< (0) ^ ^ f)}) the functions h^ and hf belong 

to the space C^([—|]), take positive values on the segment [—|, |] and h^(—|) < hi, 
h^{{) < ^ 2 - With the help of functions h^ and h^ it is possible to describe both a local 
narrowing and a local widening. 


In the domain we consider the following mixed boundary-value problem: 

( 



-Aus 

{x,y) = 

fix, 

f), (a;,|/)efi£. 


-dyUe{x 

^y)\y: 

=±eA = 

h) 

[x), X e Ie \ i = 1,2, 



Uei- 

-1,2/) = 

0, 

2/ e 

(1) 


Us 

(1,2/) = 

0, 

2/ e Tf\ 



duUe 

{x,y) = 

0, 

ix,y) G Fe, 



= ( 

1), 

dy is 

i the outward normal derivative, and the 


boundary of the joint is described by the formula 


n = asJA U (F® X ® F’)) 


2=1 


Assume that the given functions / and are smooth in the corresponding domains 

of dehnition. 

It follows from the theory of linear boundary-value problems that, for any hxed value of 
e, problem ([1]) possesses a unique weak solution from the Sobolev space such 

that its traces on the vertical end sides of the domain fig are equal to zero, i.e., Ue\x=±i = 0, 
and the solution satisfies the integral identity 

2 


Vug • Vf) dx dy = 


f^pdxdyT£Yl /(i) 

i=l 


V’ dx 


( 2 ) 


for any function -0 G such that 'il>\^^j-i = 0. 
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Remark 1.2. In the right-hand side of identity we introduce the abridged notation 


E L V' * -E E j II i>dx + Bj2 I 

i=l '' i=l ic j=l -J is ' 


^|J dx 


and use it in what follows. 


The aim of the present paper is to construct and justify the asymptotic expansion of the 


solution Me as e —)• 0 . 


2 Formal asymptotic expansion 

2.1 Regular part of the asymptotics 

We seek the regular part of the asymptotics in the form 


M® ■ = 
“oo • 


+ 00 

, {x,y)eQ^f\ i = l,2. 

k=2 


( 3 ) 


Formally substituting the series (|3]) into the differential equation and into the hrst boundary 
condition of problem ([1]), we obtain: 


+ 00 


- ^ h)-Yl 


+00 


k-2f)2 ,,(*)/ 


+cx> 


k=2 


k=2 
+ 00 


k=2 


J2e^dr,u^f\x,r]) 


k=2 


^2, ,h) 

.k-2 ^ ^k 

da;2 


2 (i) f \ 

e 


y 

'x)^f{x,r]), V = -^ 


'n=±- 


where dx = dldx, = d'^fdx"^, drf = dldr], 8“^.^ = 8“^/drf. 

Equating the coefficients of the same powers of £, we deduce recurrent relations of the 
boundary-value problems for the determination of the expansion coefficients in ([3]). Let us 
consider the problem for u ^2 • 


(hr™ 


-8^uf’(x,T]) = f{x,7]) + 


-8,^uf{x,ri)\^^ = (p±(a:), xG# 

II zsz 2 

{u^f\x,-))ri = 0, xe#. 


(x), T] G Ti, 


( 4 ) 


Here T* = , (m(x, -jjxi := m(x, r])dr], i = 1,2. 

For each value of i, the problem (jl]) is the Neumann problem for the ordinary differential 
equation with respect to the variable r; G T,; here, the variable x is regarded as a parameter. 
We now write the necessary and sufficient conditions for the solvability of problem (jT]) and 
obtain the following differential equation for the function uj^'^ : 


hi 


d^ujf 

dx^ 


(x) = / f{x,rf)dr] — g)^^{x)-\-i)’f}{x), x G(i = l,2). 


(*)/ 


( 5 ) 


'Ti 
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Let be a solution of the differential equation ([5]) (boundary conditions for this differential 
equation will be determined later). Then the solution of problem (jlj) exist and the third 
relation in (jl]) supplies the uniqueness of solution. 

For determination of the coefficients i = 1,2, we obtain the following problems: 

{uf{x,-))r^ 


= 0, xe# 

= 0 , xe 


Repeating the previous reasoning, we find = 0 and —r^(^) = 0 a; G /i*\ i = 1,2. 

dx^ 


Let us consider boundary-value problems for the functions Uu\ k > 4, i = 1,2 : 


-d^^u^j^\x,r]) = ^^-^{x) + dl,^uil2{x,T]), ijeTi, 


dvut\x,r])\ = 0 , xel 


(i) 


'ri=± 


(7) 


{u^\x,-))Ti = 0, xel^"’. 


(*) 


Assume that all coefficients ^ • • • > • • • > of expansion ([3]) are determined. 

We find and from problem ([7]) . It follows from the solvability condition of problem 

([7]) that 


dx^ 


dluk-2{x,r])d7] = -dl ull^{x,vi)dri^ = 0 , 


i.e., is a linear function solving the differential equation 


d^uf 

dx'^ 


{x) = 0 , 


X e 


( 8 ) 


Remark 2.1. Boundary conditions for the differential equations ^ and I® are unknown 
in advance. They will be determined in the process of construction of the asymptotics. 


Thus, the solution of problem ([7]) is uniquely determined. Hence, the recursive procedure 
for the determination of the coefficients of series ([3]) is uniquely solvable. 

Remark 2.2. By using the recursive procedure for the boundary-value problem one can 
easily show that the functions identically equal to zero for odd fc = 2p-|- 1, p G N. 


2.2 Boundary asymptotics near the vertical sides of domain 

In the previous section, we have considered the regular asymptotics taking into account the 
inhomogeneity of the right-hand side of the differential equation in ([T]) and the boundary 
conditions on the horizontal sides of the thin domain In what follows, we construct 
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the boundary part of the asymptotics compensating the residuals of the regular part of the 
asymptotics at the left side of and the right one of . 

At the left vertical part of the boundary of we seek the boundary asymptotics for 
the solution in the form 


+ 00 




( 1 ) /1 + X y 


k=0 


( 9 ) 


Substituting ([9]) into ([T]) and collecting coefficients with the same powers of £, we obtain 
the following mixed boundary-value problems: 

0 , 

0 , .. 

fl'l ('1'i 

nl"(0,r,)= ,,6 Ti, 


{^,ri) e (0,+cx)) X Ti, 
^ e ( 0 , +cxd), 


0 , 


^ -Fcx), veTi, 


where 




1 -I- X 


, V = ^, = -a;[!j2(-l), A: = 0,1, 


£ £ 

.,(1) 


^k’iv) = - 4 ^ 2 (- 1 )> k>2. 

Using the method of separation of variables, we determine the solution 


+00 


p=0 


_ 2 e 2 l ^ / 2p7r 


a^e H ’ cos 


. (2p+l)7r ^ 


(2p l)7r 
sm 1 — - - 7] 


of problem ffT0|) at a hxed index k, where 


( 11 ) 


^ 

h 


- J ^k\v) cos dr], 6^ dy, 


tiri — 


h 




It follows from the fourth condition in flTOl) that coefficient must be equal to 0. As a 
result, we arrive at the following boundary conditions for the functions {^^.+ 2 } : 


'^i+2(“l) “0) ^ ^ 1^0- 


( 12 ) 


f 

At the left vertical part of the boundary of kll , we seek the boundary asymptotics for 
the solution in the form 

r(2) 

e ’ 5 


+00 


nW:= 


k=0 


(13) 
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(<2'\ 

We obtain the following problems for the determination of the coefficients {11^ 

h) = 0, (r,??) e (0,+oo) X T 2 , 

-dr,u^;\C,v)L-±!^= 0, re(0,+oo), 

ni‘’(0,i()= «?>(.;), .jeT.,, 


0 , 


—)■ + 00 , T] G T 2 , 


( 14 ) 


where 


r = h = |, fc = 0 , 1 , 


^k\v) = -4?2(1)> k>2. 

Similarly we find the following solution of problem (IT^ at a fixed index k : 

+ CXD 


nf «*.<;) = E 

p=0 


(2\ _2E2Lf* 

a^p ‘e '*2 cos 


2p7r 

^2 


(2p+i)^ c* . /(2p+l)7r 


r] j + '‘2 ’ sm ( ' ~ ^ ^— r] 


where 


= — 
^ h2 


j ^k\v) cos dr], ^ 


^2 

2 


(15) 


^ j 't’t'(l) sin ) *5 

h.9 


2 


2 


2 


.( 2 ), 


.( 2 ) 


.( 2 ) 


^0^ = ^ / ^ y -^ 1 + 2 ( 1 ) = - 4 + 2 ( 1 )- 


It follows from the fourth condition in flT4)) that the coefficient is equal to 0. This is 
possible if 

4 + 2 ( 1 ) = 0, fee Nq. (16) 


4b(i) = 

I - O -f/^rr^ 1/^ - Oon I T fY\ /o 7/> I/V/o 


Remark 2.3. Since u\ = 0 for k = 2p + l^ p , we conclude that = 0 and, hence 

n? = 0, = 0) pen, i = i,2. 

Moreover, from representation m and l(T^) it follows the following asymptotic relations 

ni'lf.’l) = 0(exp (—as ( ^+«>, 

nf(e.>)) = 0(exp(-^r)) as r^+«s. 


(17) 


Equalities flT^ and flT6l) specify the boundary conditions at points —1 and 1 for all 
functions and }, respectively. 
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2.3 Inner boundary part of the asymptotics 

To obtain conditions for the fnnctions and at the point 0, we introdnce an 

additional internal asymptotics in a neighbonrhood of the joint. For this we pass to the 
following variables ^ = f and rj = K Then forwarding the parameter £ to 0, we see that 
the domain is transformed into the unbonnded domain S , which is the nnion of joint 
and two half strips = (—oo, —|) x Ti, = (|,+oo) x T 2 , i.e., S is the interior of 

s(b u “(0) u s( 2 ) . 

Let us introduce the following notation for parts of the boundary of the domain S: 

• is the horizontal parts of the boundary i = 1,2, 

• T = dE\ U dE2'>) . 

We seek the inner expansion in the form 

+ OD 

= ( 18 ) 

fc=i ^ ^ 


Substituting (ITSll into (II]) and equating coefficients at the same powers of £, we derive the 
following relations for {W} : 


where 




(e,h) e s. 


duNk{^,7]) = 

0 

K.>))6r, 


k{^,V)\r,=±!^ = 

dho. 

(-i)*e e ( 

|,+oo), i = l,2. 

Nk{^,v) - 

'^ 1 + 2 ( 0 ) + '^k 


+ 00 , r] eTi, i = 

III 

III 

JZ5 

Fk{^,v) = 

ck—2 ^^—2 -f 

(A:-2)! 9x^-2^°’"^^’ 

(e,h) e S, 


(19) 


SW ^ ^ 0, 


(-l)*ee(|,+oo), z = l,2. 


^k-2 


vl/«(^,,^) = ^^(0), ^ = 1,2, 




( 0 ) + ^-r^( 0 )+E * = 1,2, k>2. 




( 20 ) 


dx 


k\ dx^ 


j=o j\ dx^ 


The right hand side and boundary conditions for problem (IT^ are obtained with the help 
of the Taylor decomposition of the functions / and at the point x = 0. The fourth 
condition in flT^ appears by matching the regular and inner asymptotics in a neighborhood of 
the joint, namely the asymptotics of the terms {W} as ±cxd have to coincide with the 
corresponding asymptotics of terms of the regular expansions (I3|) as x —)■ ±0, respectively. 
Expanding each term of the regular asymptotics in the Taylor series at the point x = 0 and 
collecting the coefficients of the same powers of e with regard to (0), we get relations (EOj). 
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A solution of problem flT^ is sought in the form 

Nk{^,v) = ^k\^,v)xi{0 + ^k\^,v)X2{0 + Nk{^,v), 

where Xi ^ 0 < Xi < 1 and 


Xi(0 = 


0, if (-i)*e < 1 + |, 
1, if (-l)*e>2 + |, 


i = 1,2. 


Then has to be a solution of the following problem: 

-A^r,Nk{^,r]) = Fk{^,T]), (^,r/)eS; 

dM^,v) = 0 (e,r/)er; 

— R(b 


(-!)■{ e(^,+a)), s = l,2, 


Wid(f,i)) -> TblO). (-1)’^-»+0O, i; 6 Ti, i = l,2, 


where Fq = 0, 




du. 


(b 


Fk{^,v) = Y1 


2 = 1 


(« 


2=1 

,(b 


dx 

k-2 ai„.d) 


d(,X'> 

®^fc+l2n^ I ^ ® '^2 I W "(c\ 

^ ^ d)m (0 


dx 


J=0 


j! 9x1 


+2 


' ‘"(0) + + E ^^(0.,) )x'K) 


V dx 


(fc — 1)! dx^ 


ck—2 f)k—2 f 

^ ^(o,d)x.(0 


(fc — 2)! dx^~‘^ 


^ (j - 1)! 9x1 

tfc —2 ^^“2 £ 

^ (A:-2)! 9x^-2^°’"^^ 


and 


tA;-2 Jk-2,M) 

?(b_5(b_n '5 « ‘di 


( 21 ) 


( 22 ) 


BJV = Bi‘( = 0, 


BEK) = 


(fc —2)! dx^“2 


^( 0 )(l-x*(O), * = 1,2, A) >2. 


To study the solvability of problem ([22]), we use the approach proposed in [8]. Let 
be a space of functions inhnitely differentiable in S and hnite with respect to i.e., 

VneCo“^(S) 3F>0 V(e,r/)eS \^\ > R: v{^,v) = 0- 


o“( 


We now dehne a space R ;= (C^(S), || • H-t^) , where 


l^^lk = '\/ / |Vn(^,r7)|2d^dp+ / |n(^,d)Hp(0N^c(d, 


and the function p(^) = (1 + |^|) ^ G M. 
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Definition 2.1. A function Nk from the space H is called a weak solution of problem 
if the identity 


j ViVfc • Vvdfdr] 


+ 00 


Ftvdidr,^ I 


;( 2 ), 


( 23 ) 


—oo 


1 

2 


holds for all V . 


From lemma 4.1, remarks 4.1 and 4.2, corollary 4.1 (see 0 ) it follows the following 
propositions. 

Proposition 2.1. Let p~^Fk G F^(S), G L^(|,+oo) and G L^(— oo,—|). 

Then there exist a weak solution of problem (d^ if and only if 

^ n d(;dri T ’ b£>( 0 T Bg(?) if = 0. (24) 


This solution is defined up to an additive constant. The additive constant can be chosen to 
guarantee the existence and uniqueness of a weak solution of problem (d^ with differentiable 
asymptotics 


Nk{^,v) 


C’(exp(^ 0 ) as -oo, 

+ C’(exp(-^ 0 ) as +00. 


(25) 


Proposition 2.2. The corresponding homogeneous problem for problem ^2E) 


— = 0 in E, = 0 on dE, 


(26) 


has a solution OIq that does not belong to the space Ti and it has the following differentiable 
asymptotics: 


%(e,h) 


+C>(exp(^ 0 ) as ^ ^ -oo, 

<^0 + +C>(exp(-^0) as ^ ^+00. 


(27) 


Any other solution to the homogeneous problem, which has polynomial growth at infinity, can 
be presented as a linear combination ai + ooffio- 


Proposition 2.3. If the domain E is symmetric about the horizontal axis, the function 
is even with respect to the variable p {Fk is odd with respect to p) and = —B^^, i = 1,2 

i^kl = ^kli * ~ I 52 ), then solution Nk is an even (odd) function with respect to p. If Nk 
is an odd function, then the constant in (d^) is equal to zero. 


Remark 2.4. Using the second Green-Ostrogradsky formula, similarly as was done in re¬ 
mark 4 ■ 3 (W, constant S)) {k G N) in (d^) can be found as follows 



+ 00 


fFk{^,p)dfdpT I ^B^kliOd^T I I dfif) Nk{^,v) da^. 


( 2 ), 


( 28 ) 
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It follows from Proposition 12.21 that problem (12^ at A; = 0 has a solution if and only if 


4i)(0)=a;r(0); 


(2), 


(29) 


in this case 


N„ = N„ = w<‘>(0). 


(30) 

Let us verify the solvability condition (|2^ . Taking into account the third relation in 
problems (jl]) and ([7]), the equality fl2TD can be re-written as follows: 


h- 


du 


( 2 ) 


du 


( 1 ) 


dx dx 


( 0 ) 




-\-hi 


{k — 1)! dx^ 


(0) Xi(Oc?^ + ^2 




(A; — 1)! dx^ 


( 0 ) X2{Od^ 


_ 1_9 
2 ^ 


+1 


hi 




+ 00 


(A: —2)! dx’^ 


(0)(i-xi(O) - h 






{k — 2)\ dx^ 


■(0)(l-X2(O) di 


ck—2 f 

r^{0,T]) d^dr] = 0, A: G M, k>2. 


y=(o) (A: — 2)! dx^~‘^ 

Whence, integrating by parts in the hrst two integrals with regard to (jS]), we obtain the 
following relations for {u^^} : 


du 


( 2 ) 




dx 


i0)-h 


du 


( 1 ) 


^ dx 


-(0) = -dl, k eN, A: > 2, 


(31) 


where d 2 = 0, 




2=1 


(.-2)1 


+ 


^ didii, ken, k>3. (32) 


k{o) {k — 3) \ dx^ ^ 


Hence, if the functions and uf'* satisfy ([SI]), then there exist a weak solution of the 
problem fl22|) . According to Proposition 12.11 it can be chosen in a unique way to guarantee 
the asymptotics (|25|) . However, we do not take into account the limit relations at inhnity 
in fl2^ (see the forth condition). In order to satisfy them we add n;^^2(0) to our solution 
(Proposition 12.11 gives us that possibility) and derive the following conditions: 


d‘’(0) + '5t2 = T’(0) . ‘'eN, .>3. 


,( 2 ) 


(33) 
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As a result, we get the solution of the problem (IT^ with the following asymptotics: 


+ C’(exp(^0) as ^ ^ -cx), 




Nk{^,v) = 


Let us denote by 


d+2(0) +'■'t'K.'))+ O(exp(-^0) as ^-»+oo. 


,< 2 ) 


( 34 ) 


,,, *:eN. 

4+2(0) +h’K.’rt. ^>0, 

Remark 2.5. Due to functions {Nk — Gk}keN o,re exponentially decrease as f ^ ±cxo. 




2.4 Limit problem 


Relations fl2^ . fl33|l together with (jST]), (|5]), (jH]), flT^ and ffT6|) complete boundary-value 
problems to determine the functions 

So for the functions and that form the main term of the regular asymptotic 
expansion (|3]), we obtain the following problem: 


n2 (b ^ 

-hi I (x) = Fhl{x), X e A, i = l,2, 
dx^ 

= n;f(0), 


duj. 


( 1 ) 


duj. 


( 2 ) 


(35) 


ft,=;^(o) = 


dx 

.(!)/ 


dx 


= 0, uf\l) = 0, 

where h = (-1,0), h = (0,1), 


Fh\x)-.= / f{x,7])d'r] — (p^jf{x)F‘p_{x)^ x G A, i = l,2. 

Jr, 


(36) 


The problem 0351) is called limit problem for problem ([T]). The solution to 0351) is given by 
the following formulas: 


CJ. 


( 1 ) 


(x) = — / (s — x)F‘^^\s)ds 


hi 


-1 


(x -b 1) 
hi + ^2 


n {^^s-l)F^^\s)ds+ f{l-s)F^^\s)ds'), xe/i; (37) 


-1 
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u. 


( 2 ) 


(x) = — / (s — x)F‘^^\s)ds 


h 


1 u 


For next functions : A; > 3}, the problems take the form 


-h, 




x) = 0, 


dx'^ 

^0) = ^k\^) ~ ^k-2^ 


X G li, i = I, 2, 


du 


( 1 ) 


du 


( 2 ) 


(39) 


Ai^(O) = A,^(0) + <iJ, 

dx dx 

k‘’(-i) = 0. k"’(i) = 0. 

It is easy to verify that the solution to problem fl5^ is given by the formulas 

= - h~Ph - (^+ 1 )’ 


(40) 


3 Complete asymptotic expansion and its justification 

From the limit problem fl35|) we uniquely determine the hrst term of the asymptotics Ci ;2 of 
series Q- Next from the equality fl3Up we obtain the hrst term Nq of the inner asymptotic 
expansion fITS]) . Then we rewrite problems (jl]) in the form 


-dl^uf{x,ri) = f{x,7]) - h. ^FW(x), t] e T^, 

ii)l 


^ = (p^^\x), X e li 

2 


* = 1,2, 


(41) 


(** 2 ^, •))t = 0 , X e li, 


and hnd that 


U 2 \x,r]) = — / ^ {r] — t)yf{x,t) — h^^F^'^\x)jdt — r](p^l\x) + a 2 \x), 
where function are uniquely determined from third condition in fjTTll . i.e. 


( 42 ) 


a? (a;) = 



{rj — t) f{x,t) dt dr] — Q ^h‘^F^^\x), i = l,2; 


/ 'Y' / '^1 

i-z 2 
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functions and are defined by relations (l36D . 

Now with the help of formulas ffTTj) and flT^ . we determine the hrst terms and 
of the boundary-asymptotic expansions ([9]) and flT^ respectively, as solutions of problems 
flTOl) and flT^ that can be rewritten as follows: 


= 0, (^,77) e (0,+cx)) X Ti, 

-drj^?{^,v)\r,=±^ =0’ (0,+oo), 

>1 -i- 2 

0, ^^+CX), 

= 0, (r,h) e (0,+cx)) X T2, 

(e,h)l,=±i^ = 0, re(0,+cx)), 

4')(0,r/)= heT2, 


(43) 


(44) 


n^'^(r,h)^ 0, 


C +00, T] eT 


2- 


The second term Ni of the inner asymptotic expansion flTS)) is the unique solution of the 
problem (IT^ that can now be rewritten in the form 


-A^riNi{^,r]) = 0, 

d,N,{^,v) = 0, 

d*^ + 


Ni{^,v) 


hi -|- h2 




(e,h) e s, 

(C, //) e as, 

- 1 )*^ -f-cx), 77 e Ti, 7 = 1,2, 


with asymptotics fl34ll . Recall that the constant dg is determined by formula 
constant Sf is also uniquely determined (see Remark I2.4p by formula 


(45) 
and the 


6t = —^ Ni{^,r]) da^n^ . (46) 

Thus we have uniquely determined the first terms of the asymptotic expansions ([3]), (jO]), 
(m and ([H]). 

Assume that we have determined coefficients , ^ 2 n- 2 i "^2 ^ ■ ■ ■ 1 '^ 2 n -2 of 

series ([3]), coefficients ..., n 22_2 of the series ([9]) and flT^ respectively, coefficients 

Ai,..., N 2 n -3 of the series flTSl) and constants hj*",..., S^ns ■ 

Then, using formulas (HOjl . we write the solution 022^-1 of problem (|3^ with the constant 
.3 in the first transmission condition. It should be noted that constants {d^}A :>3 depend 
only on / and 7 = 1,2 and they are uniquely defined by formulas ([32]). Further 

we hnd the coefficient N 2 n -2 of the inner asymptotic expansion ffTSj) . which is the unique 
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solution of the problem flT^ that can now be rewritten in the form 


/ 




\ 


-A^riN2n-2{^,V) = 

d,N2n-2{^,V) = 

-dr,N2n-2{^,V)\rj-±!li = 

'1^2 

N2n-2 ~ 


t2n—4 a2n—4 £ 

- -zi(0,h), (e,h)es, 


(2n — 4)! dx‘^^ 

0 , 

^2n-4 ^2n-4^« 


(e,h)er, 

(0), (-l)*ee (|,+oo), z = l,2, 


(2n —4)! 

d*2n-2 + {-^yh^-^dtn-A , (-1)'+'^/^-! + <^2^3 , 

H-;——;-4 


hi + /i2 

‘2n-2 j2n-2, ,(*) 


hi + /l2 


^2n-2 ^2n-2^yv ^ 5'M2n-2-i 


^-(0) + V - 

(2n-2)! 2 ^ ^ J 


{o,v), 


dx^ 

1)*^ ^ +00, rj eTi, i = 


1 , 2 , 

(47) 


and N 2 n -2 has asymptotics fl341) . 

Knowing (see (|28|1 ) and using relations (IdOil . we get the solution uj 2 n of problem 

(|3^ . Next coefficient A^ 2 n-i of the inner asymptotic expansion (1T8|1 is dehned as the unique 
solution to problem (IT^ that can be rewritten in the form 


-A^r^N2n-i{^, v) ^2n - 3 )! dx^^-^ 

duN2n-l{,^,V) = 0, 

^2n-3 ^2n-3^W 


c2n—3 pj2n—3 f 

^ r^(o,h), (e,h)es. 


dr)N2n—li^-i V) L=4-li 
>1-1-2 

N2n-l ~ 


( 0 ), 


(e,h)er, 

(-i)*ee (i+oo), z = l,2. 


(2n —3)! (ia;2"'“3 ^ - V 2 

dln+l + {-^yh^-^dtn-l , (-ir+'dL + /^ 3-4 52V2 


hi + ^2 

J2n-1, ,(*) 


+ 


hi + ^2 




^2n-l ^2n-l^W ^ 2^3 


(2n — 1)! (ix2" 


^(o)+ E 7 


dx^ 


{o,v), 


Coefficients U 2 I, i = 1,2, are determined as solutions of the following problems: 


j=o r- 

(- 1 )*^ -t +00, 7 e Ti, i = l, 2 , 

(48) 


-diyUix, r]) 

= 

;3 .!!: 
1 

to 

£(a:,h)l^=±K 

= 0, 

X e li. 

(4n(^> •))t, 

= 0, 

X e h. 


(49) 


We note that solvability condition for problems fl49|) takes place, because (M 2 n- 2(^5 ■))Ti = 0, 

i = 1,2. 

Finally, we hnd the coefficients and fl^^^ of the boundary asymptotic expansions 
([9]) and flTSD respectively as solutions of problems flTU]) and ffT41) that can be rewritten in the 
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form 

(^, ??) = 0, (^,? 7 ) e (0,+cx)) X Ti, 

-^r,n^n(^,r/)| ^ =0, ^ e (0,+Oo), 

I fl—nz 2 

nS(0,?/)= -MS(-l,r/), r/eTi, 

. ^ 0) ^ +00, r] G Ti, 

-Ag*,n£(e,^) = 0, {C,v)e{0,+oo)xT2, 

-a^ng(e,r 7 )|_^ = 0, re(0,+cx)), 

/ '1^2 

n?j(0,??)= -42(1,r?), T^eTa, 

. ^•2ni^*^v)^ 0) +00, r; G T 2 . 

Thus we successively determine all coefficients of series (Ej), (Ej), flT^ and flTSll . 


( 50 ) 


(51) 


3.1 Justification 

Let us introduce the following notations 


^ y\ 

Uk{x, -) = 


,d) 


“r 44) ’ ^ < 0 


ujI^\x), x<0 
,( 2 )/ 


u 


( 2 ) 


, (Xk{x) = 


44)’ ^>0 

and define the coefficients of regular asymptotics as follows: 

=(= 

Ui 


ui^ix), a; > 0 


A: G N, k>2 


= Uk + a;fc+2(x), A: G No (^^o = mi = 0). 

With the help of the series (E]), ([9]), flT^ . flTSl) we construct the following series 



+ 00 

+ 5 :'’“ 

k=l 




+ x+(x)n 


(2)4 


X 


2k 



{x,y)ene, (52) 


where a is a hxed number from the interval (|, 1), Xh smooth cut-off functions 

defined by formulas 


Xi{x) 


1, if |x| < /, , 

0 . if|x|> 2 i. 


1, if |1 T 3;| < 5, 
0, if |1 =F x| > 2(5, 


and 5 is a sufficiently small fixed positive number. 
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Theorem 3.1. Series fl 52 p is the asymptotic expansion for the solution of the boundary- 
value problem ID in the Sobolev space i-e., 

VmeN 3Cm>0 3eo>0 V£e(0,£o): he - (53) 

where 

2m 

utHx,v) = J2^'‘ 

k=0 
m 

k=l 

is the partial sum of fl 52 p . 

Remark 3.1. Hereinafter, all constants in inequalities are independent of the parameter e. 

Proof. Take an arbitrary m G N. Substituting the partial sum in the equations and 

the boundary conditions of problem ([1]) and taking into account relations fl3^ - fl5TD for the 
coefficients of series ([52]), we hnd 



where 


D 

+ / (x, f) =^R^>(x,p) =: R<r\x,y). 



^2771 




d^U2.,r / y\ 

dx^ \’e)’ 


(55) 


(56) 


d”h.9) = - dEGtii.r,)) 


k=l 


+ £ 


-2afXl , 
dC 




(57) 


(:=— B=- 71=^ 


= (2^ '^(h^«n(^fc(e,r7) + ^^(a;)n^l^(e,r7) 


d^X~ 


k=l 


2e 


dx 
-idx 


dx 




X 


e°‘J (2m — 2)! 


dx‘^ 

ix)d^^2k (^> V) + (^> V) 

X - Q2m-lf 






, (58) 


7 (^’ f) 
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(2m - 1)! 


X — z 


2m—1 


(P^+^Uo 


— 1 ^{l-a) 2 j ^ _ .2m-2j-l g 2 m- 2 j+l 


-E 


^ (2m - 2j - 1)! 


^^2m+l 

“2i / y 


z') dz 


Q^2m—2j-\-l 


[z, I) dz 


2f(0k.^. 


(60) 


„,k pO(2m—1) j ^(1—a)2m—a^ *^^2m+2 




(ix 


1 


(2m)! 


2m 


^2m+] 


022 


^^2m+l 


z) (iz 


E 


M-Q)2j r 

^(2m-2j)! J 
From fl56l) we conclude that 




3Cm>0 3 £o>0 V£G(0,£o): sup R^^i{x^y) 

(x,y)en, 


< (62) 


Due to the exponential decreasing of functions {Nj. — (see Remark 12.51 and 

(1T71F and the fact that the support of the derivatives of cut-off function xi belongs to the 
set {x : le^ < |a:| < 2Ze:"}, we arrive that 


sup 

{x,y)&Q.e 

similarly we obtain that 




< CmS 


— 1 — Q. 


ttI 


max (hi, ^ 2 ) e^~ 


sup 

{x,y)ene 


R‘fl\x,y) <Cme 7xp 


7i6 


max (hi, h 2 ) s 


(63) 


(64) 


We calculate terms R^J^\ j = 4, 5, 6 with the help of the Taylor formula with the integral 
remaining term for functions / , 022 and {M 2 fc} at the point x = 0 . It is easy to check that 


sup 

{x,y)Gne 




j=4,5,6. 


( 65 ) 


The partial sum leaves the following residuals in the boundary conditions: 

dyU^'^\x,y)\^^^^^+eip^^\x) = R^™^^(x), x e #, * = 1,2, 

U^^\±l,y) = 0 , yeT^\ * = 1 , 2 , 

dM^"^\x,y) = 0 , {x,y)eTe, 
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where 




Xi 


X 


£"/ (2m — 2)! 


X 2m-2 ,2m-2,Ai) 

^ ^ ^{Z)dz, ^ = l,2. 




It follows from fl66|) that there exist positive constants Cm and Eq snch that 


V£e(0,£o): snp R^J2l^{x) 


x£l, 


< CmE 


l-\-a{2m—l) 


, t = l,2. 


Using estimates (|62il - (|65ll and (|67|l we obtain the following estimates: 


dM 

-^£,1 


L^fle) 


< Cm \/hl + ^2 2 ^ 


R 


(m) 


■e,2 


L2(n^) 


R 


(m) 

'£,3 


< Cm\/l max (hi, ^ 2 ) £ 2 exp 

< Cm\/l max(hi,h 2 ) 5^e~^ exp ^ 


tt/ 


max (hi, ^ 2 ) 

TT^ \ 

max (hi, h 2 ) £/ ’ 


R 


,(m) 


'£,4 


L2(0,) 


< ( ihi + llh2 + 


R 


(m) 




L^iQe) 


<CmVlh+lh2 J=5,6, 




,{m) 

'e,(*)± 




<„, e'+“P"-i>, * = 1,2. 


( 66 ) 

(67) 

( 68 ) 

(69) 

(70) 

(71) 

(72) 

(73) 


(74) 


Thus, the difference Ws '.= — f/i™^ satishes the following system: 

-AlU, = infi„ 

-a,fU,(a;,±£f) = R^^)^{x), a:e#,* = l,2, 

1U,(±1,2/) = 0, |/eT«, z = l,2, 

^,,144 = 0, on Tg, 

This means that the constructed series 0521) is a formal asymptotic solution of problem (|T|). 
From fl7^ we derive the following integral relation: 

2 

f \\/We\‘^dxdy = j W, dxdy t ^ / '^£,(b± ^^l?;=±£U dx. 


i=l 


r(0 


In view of the Friedrichs inequality and estimates fl68|l - fl73|) . this yields the following 
inequality: 

/ \VW,fdxdy < Cm + C,, i) ^ 

d j=i 
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This, in turn, means the asymptotic estimate fl5^ and proves the theorem. □ 

Corollary 3.1. The difference between the solution of problem m and the solution U 2 
of the limit problem fSR) admits the following asymptotic estimate: 


\\us — uj2\\m{ne) ^ CqE. (75) 

In thin rectangles := x i = 1,2, the following estimates hold: 

3 

\\ue - UJ2 - < Cie2, i = l,2; (76) 

in addition, 

<C2£, * = 1,2, (77) 

m.Q;yi \E'f\ue){x) — 0J2{x)\ < C^e, i = l,2, (78) 

-rii) 

where /i,a := (—1,—2/e"), := (2/e",l), a is a fixed number from the interval (|,1), 

UJ 2 , is defined by the formula fl4U p and 

Ei'\ue){x) = ^ [ ufix, y)dy, i = l,2. 
s hi Jrfi 

In the neighbourhood 12® := 12^ n {{x,y) : x G {—el, el)} of the joint, we get estimates 

||Va;yM£ - Vg,, A^i||^2(f^(o)) < ||m£- a;2(0)-eiVill^i^^(o)^ <C'4ei"+i (79) 

Proof. Denote by Xifi') •= X/(^)- Using the smoothness of the functions {uJk}t =2 and the 

exponential decay of the functions {iV^ — Gfc}fc=i ^2 , U® and If® at inhnity, we deduce the 
inequality fl75|) from estimate fl5^ at m = 1: 


< 


\\us - W 2 II HHQe) - ll“= “ Whhu,) + “ ^2 + xU No + e{{l- xU) ^3 + xU ^0 

~ ^“2 + ^a) + x\,a^2 + x n® + x’^n® j ^ 

< Cie5«+2 + ||a;2 - W2(0)11^1(^(0)^ +e ||iVi||^i^^(o)^ + ll^2||j;^i(f2(o)) 

(^llx n2 IIj^i(q(i)) + llx 172 llj^i(Q(2))j 

2 

||xf,a('X'2(0)-a;2)+e^ ((1 - XIq)‘^3 + xla^l) ((1 “ xla)i'^2 + ^a) + Xt,a^2) 

i=l 

Cl e2"+2 +5Ci +e ||2 Vi||j^i(2(o)) + e^ II^2 ||hi(h(o)) (IIx U^ + llx'^U^ ^|| 


Hpnfi) 
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I 

\\xla \\XU (^2 - G'2)||j:^i(f2(i)) 

2=1 


E ll , / d(jj2 /„N d^ojo /„x \ 

||x,„ (a;2(0) + X—(0) + y ^(0) - u:^) 

2=1 


d‘^U2 




\xd,a ( 1 ^ 3 ( 0 ) + 

2=1 




Z 

(“ 2 ( 0 > •) - “ 2 ) + ||xf,a (^ 4 ( 0 ) - U^) Li(^W)) 

2=1 


+ ll^3||j:^i(t2W) + \\U2 + ^4\\H^ni^))J < C'l 

2=1 

Again with the help of estimate (j53|l at m = 1, we deduce 

\\Ue-UJ2 - ^ he “ UhllH^Q,) + e‘^\\u2 + UJ4 + X ~) 

< Cieh+h +C2e'^, 

whence we get (j76ll . Using the Cauchy-Buniakovskii-Schwarz inequality and (1761) . we obtain 

inequalities d77|) . Since the space H^{Ie%) continuously embedded in from d77|) it 

follows inequalities fl78|) . 

From inequalities 

ll'i^e “ 1 ^ 2 ( 0 ) — £ -^i|Ij 71 (q( 0 )) — he ~ Uh\\m(ne) + < Ci £2"+2 J\- 8^ 


it follows more better energetic estimate fl7^ in a neighbourhood of the joint hi 


( 0 ) 


□ 


Remark 3.2. If (p^f} = 0 and the function f depends only on the variable x, then all 
coeficient {u 2 k}-, } and are equal to 0. In this case the asymptotic series fl 52 1) 

has the following form: 


Y^h({l-Xi (y))a;fc+ 2 ( 2 :) + Xz (y) (y |) V h,y) e 
k=o \ / 


(80) 


and the residual terms {R^hj^}hi also simplified respectively, but the asymptotic esti¬ 
mates (1 54 p remain the same. Nevertheless, as follows from the proof of Corollary 1,9. il the 
asymptotic estimates (1 76 p — (1 78 p become better: 


max 

-=(i) 


\\Ue -UJ2- < Cl £2“+2 , i = 1, 2; 

\\EhiUe) -UJ2 - < C 2 i = 1, 2; 

Eh{us)ix) -U2(x) -e:w3(x)| < i = 1,2. 


(81) 

(82) 

(83) 
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4 Conclusions 


1. The energetic estimate (j75jl partly confirms the first formal result of [4] (see p. 296) that 
the local geometrical irregularity of the analyzed structure does not significantly affect on 
the global-level properties of the framework, which are described by the limit problem (|35il 
and its solution uj 2 (the leading term of the asymptotics). 

But now, due to estimates fl76|) and flHTj) - (1831) it became possible to identify the impact 
of the geometric irregularity and material characteristics of the joint on the global level (the 
second term of the regular asymptotics ([2]) depends on the constant hj*" that takes into 
account all these factors (see (IT6D B. This our conclusion does not coincide with the second 
main result of [Tj (see p. 296) that “the joints of normal type manifest themselves on the 
local level only'"'. 

In addition, in |3] the authors stated that the main idea of their approach “is to use 
a local perturbation corrector of the form eN(x/e)^ with the condition that the function 
N{y) is localized near the joint ”, i.e., N{y) —)■ 0 as |y| —)■ -|-oo, and the main assumption 
of this approach is that VyN G Li{Q^) (see (14) and similar assumptions on p. 300 and 
p. 303). 

As we see the coefficients {Nj.} of the inner asymptotics ([T8|) behave as polynomials at 
inhnity and do not decrease exponentially (see fl34l) ). Therefore, they influence directly the 
terms of the regular asymptotics beginning with the second one. Thus, the main assumption 
made in [1] is not satished. This is the second our principal disparity with results of [1]. 

2. From fl75|l it follows that the gradient is equivalent to ^ in the L^-norm over 
whole junction hie as £ —)■ 0. Since ||^||£, 2 (q(o)) = as £ —)• 0, the estimate (|75|l is not 

informative in the neighbourhood Qfj of the joint 

The form of the complete asymptotic expansion (l5^ gives us possibility to improve the 
zero-order approximation of the gradient (flux) of the solution both in the main parts /i(a, 
i = 1,2, of the junction: 

^ / N duj2, . duj'i, , 

\uAx,y) ——(x) + e——(x) as e —)■ 0 
^ ^ dx dx 

considering the geometric irregularity and material characteristics of the joint (see formulas 
fl76|) . fl8T|) h and in the neighbourhood Vtf] of the joint: 

Vm£(x,|/) ~ Vg^(A^i(C,? 7 ))|^=|,^=| as e^O 

(see (1791) b Also using estimates (153|) . we can obtain more better approximation of the 
solution and its gradient with preset accuracy. 

3. The results obtained give the right, in terms of practical application, to replace the com¬ 
plex boundary-value problem ([T]) with the corresponding simpler 1 - dimensional boundary- 
value problem (l35|) with sufficient accuracy that measured by the parameter e characterizing 
the thickness and the local geometrical irregularity. In this regard, the uniform pointwise 
estimates (|78l) and (l83|) . that are very important for applied problems, also conhrm this 
conclusion. 


23 




4. The method proposed in the present paper for the construction of asymptotic expansions 
can be used for the asymptotic investigation of boundary-value problems in graph-junctions 
of thin domains (Fig. [5]), or graph-junctions of thin perforated domains with rapidly varying 
thickness. In the last case, it is necessary to add series with rapidly oscillating coefficients 
to the regular part of the asymptotics (see m)- 
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